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Abstract. Laser trapping and interfacing of laser-cooled atoms in an optical fiber 
network is an important capability for quantum information science. Following the 
pioneering work of Balykin et al. and Vetsch et ai, we propose a robust method of 
trapping single Cesium atoms with a two-color state-insensitive evanescent wave around 
a dielectric nanofiber. Specifically, we show that vector light shifts (i.e., effective 
inhomogeneous Zeeman broadening of the ground states) induced by the inherent 
ellipticity of the forward-propagating evanescent wave can be effectively canceled by a 
backward-propagating evanescent wave. Furthermore, by operating the trapping lasers 
at the magic wavelengths, we remove the differential scalar light shift between ground 
and excited states, thereby allowing for resonant driving of the optical D2 transition. 
This scheme provides a promising approach to trap and probe neutral atoms with long 
trap and coherence lifetimes with realistic experimental parameters. 
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1. Introduction 

The development of a matter-light quantum interface using cold atoms and optical 
fibers has been an active field of research over the past years [1]. Recent advances 
towards this goal include the observation of electromagnetically induced transparency 
and the loading of ultracold atoms in hollow-core optical fibers [21 El HI |5] , as well as 
the trapping and probing of atomic ensembles via the evanescent fields surrounding 
tapered nanofibers [HI [71 El |9]. While prominent examples of off- resonant interaction 
between evanescent waves and matter have used a plane dielectric geometry for atom 
optics and interferometry [TUl E] as well as for surface traps of quantum degenerate 
gases [121 1131 E], recent progress of atom-light interactions with optical waveguides 
[H El m El E] sets the stage for the fiber integration of free-space quantum systems in a 
quantum network via quantum-state transfer between matter and light [IHl [161 [13 HE] 
and for strong coupling of single atoms and photons trapped near microcavities 
[13 EHl [211 [221 [231 [21] ■ Furthermore, these effective 1-dimensional systems may be 
applied for investigating quantum many-body phenomena in low dimensions with long- 
range interactions mediated by the waveguide [2S1 1251 [2Z1 12H] • 

One major drawback of many optical traps is that spatially inhomogeneous energy 
shifts U{r) generally depend on the atomic electronic state, limiting long-lived trap and 
coherence times necessary for repeated coherent operations [29]. This is traditionally 
alleviated by constructing a state-insensitive optical trap designed to decouple atomic 
transition frequencies from the spatially varying potential of each electronic state [50] . 
Specifically, at the "magic" wavelength conditions, the differential response of the 
dynamic scalar polarizabilities a^^\uj) for the ground and excited states a^'^^a^'^^ at 
the optical frequency u can be tailored such that both levels are perturbed identically 
with = a^^y This leads to a vanishing differential atomic level shift dUgcaiav = 
[30l [311 [321 [331 . Differential shifts for the hyperfine ground states can be minimized by 
using far off-resonant beams, whereas Zeeman coherence can be conveniently protected 
by using linearly polarized light in which the vector light shifts are zero. 

Although such magic wavelengths can be used for nanofiber traps [35] , the strongly 
guiding nature of the waveguide inevitably leads to non-negligible longitudinal electric 
fields Ez in the evanescent region, which are out of phase with the transverse field 
E_L = (Ex,Ey). Here, z refers to the direction parallel to the fiber axis, while x and y 
are the coordinates perpendicular to the fiber axis. The resulting local polarization 
at location r is in general elliptical even for linearly polarized input beams, and 
induces vector shifts ?7vector- The differential vector shift 5f/vector in turn manifests 
itself as a "fictitious magnetic field", leading to inhomogeneous Zeeman broadening 
[36] . Furthermore, the spatially varying elliptical polarization of the evanescent field 
on a scale 6r < X renders it difficult to cancel 5t/vector using bias fields, resulting in 
increased heating rate [21] and limited coherence time |37j . 

Building upon the recent realization of a nanofiber trap as proposed in Ref. [5S] and 
demonstrated in Refs. [HI [39], we propose a promising strategy for a state- insensitive 
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evanescent field trap. Differential scalar shifts 5f/scaiar between l^f) and |e) are canceled 
using "magic" wavelength conditions. The inhomogeneous Zeeman broadening SU^ector 
caused by a forward propagating blue-detuned field E'^^'*'^) is canceled by a backward 
propagating field E'^'^™'^) with a small relative frequency detuning Sjb- Thus, our scheme 
can compensate for the light shifts of the strongly guided evanescent waves to the first 
order in the space external to the dielectric fiber, leading to favorable parameters for 
the realization of a long-lived fiber-integrated quantum memory and resonant coupling 
to ultra-high quality micro-cavities based on optically trapped atoms. 



2. A state-insensitive nanofiber trap 

In this section, we discuss an ab initio calculation of the optical nanofiber trap for atomic 
Cesium. We show that the light shifts caused by the elliptically polarized components 
of the fiber's evanescent field are not negligible. We then propose a scheme to cancel 
these shifts and generate a two-color, state-insensitive, three-dimensional trap for Cs 
atoms along the nanofiber. 



2.1. ac Stark shift Hamiltonian 



We start by considering the Hamiltonian for an atom interacting with an electric field 
E in the dipole approximation: 



-d-E, 



(1) 



where d is the electric dipole operator and E is the electric field operator. Taking 
into account the atomic hyperfine structure, this Hamiltonian can be decomposed 
into its spherical tensor components parameterized by the dynamic polarizability a{uj) 

mmm-. 



= Hq + Hi + H2 
= a(o)E(-) ■ E(+) 

+ia^^>^ p 



(2) 



F(2F~1) 



l(F F + F F ) - -F'^5 



where a*^"-*, a*^^-* and a*^^^ are the scalar, vector and tensor atomic dynamic polarizabilities, 
E(+) and E^") are the positive and negative frequency components of the electric field, 
F = I + J is the atomic total angular momentum operator, with I and J the nuclear 
and electronic angular momentum operators, /i, G { — 1,0, 1} are components in the 
spherical tensor basis, and Hq, Hi, and H2 are the terms associated with the scalar, 
vector, and tensor light shifts, respectively. The light shifts ^/scalar, f^vector, and f/tcnsor 
arising from each term have been expressed explicitly in Refs. jlQl 1^ . 

For two-level atoms with ground and excited states \g), |e), the scalar shift ?7scaiar 
can be approximated by t/scaiar oc |Ep/5 for detunings 6 = u — Ua large compared to the 
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excited state decay rate F, where 00 is the electric field angular frequency and uja is the 
\g) — )■ |e) transition frequency. The ground state will experience a repulsive potential 
for blue-detuned {6 > 0) electric fields, and an attractive potential for red-detuned 
(S < 0) electric fields. The scalar dynamic polarizability a*^*^^ is in general different 
for the states \g) and |e) resulting in a differential scalar shift and a mismatch of the 
ground and excited state potentials. For the typically anti-trapped excited state, near- 
resonant driving of the transition by an additional beam with frequency U2 — oJa can 
cause significant heating of a trapped atom [29]. This situation can be remedied by the 
use of "magic" wavelengths for which a\^^^ = a^'^^ [30l |3T1 [321 

The vector term Hi of Eq. ^ induces a Zeeman-like splitting proportional to a 
projection of the total atomic angular momentum F and arises from a so-called "fictitious 
magnetic field" proportional to the ellipticity of the electric field In the case of a 
free-space plane wave propagating along the z axis. Hi can be expressed in terms of the 
Stokes operators S = (5*0, 5^., Sy, Sz) as [3T] : 

^iOC««Hey, (3) 

where e = {Sz)/{So) = j^^^piiljf "jp is the ellipticity of the electric field. For an 
elliptically polarized beam, the vector shift can be as large as the scalar shift, and can, 
for example, be used to cancel the differential light shifts of Rubidium atoms confined 
in a 3D optical lattice |13] . 

The last term H2 in Eq. ^ represents the tensor shift. It vanishes for atoms with 
total angular momentum F = 1/2 [H]. In the case of the D2 transition of Cs, that we 
consider here, it will depend only on the electronic angular momentum J for detunings 
large compared to the 6P3/2 excited state hyperfine structure, and vanish for J = ^ 
[ini SI] ■ It will therefore only act on the excited state of the Cs D2 transition, inducing 
shifts on the Zeeman rriF' sublevels proportional to m|,,. 

2.2. Evanescent optical traps using the fundamental mode of the waveguide 

When the radius a of an optical fiber is reduced well below the propagating field 
wavelength A, the resulting cladding-to-air waveguide supports only the "hybrid" 
fundamental mode HEn [71 HI]. In this strongly guiding regime, a significant fraction 
of energy of the HEn mode is carried in the form of an evanescent wave outside 
of the nanofiber. The evanescent field intensity is azimuthally asymmetric when 
the input polarization is linear [71 HI]. Fig. [l^ shows the electric field intensity 
lE'P = l-ExP + l-E'j/P + |-^2p in a plane transverse to the fiber for a single, linearly- 
polarized beam. The unit vectors [ex,ey,ez) form the basis of the {x,y,z) frame, and 
(r, 0) are the cylindrical coordinates in the transverse plane {x,y). 

By appropriately combining blue-detuned and red-detuned fields Ered and Ebiue 
in an optical nanofiber, an atomic trapping potential can be engineered from the 
evanescent electric fields [38]. Different configurations can be used for obtaining 3D 
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confinement. Here, we consider the schemes illustrated in Fig. [TJ for an infinite Si02 
cylindrical waveguide of radius a = 250 nm. In all configurations, the beams are linearly 
polarized at the waveguide input to ensure azimuthal confinement for trapped atoms. 
A pair of x-polarized red-detuned beams generates a ID lattice along the fiber axis for 
longitudinal confinement. Fig. shows the configuration used in Ref. p|, namely a 
pair of red-detuned, x-polarized beams, and a single blue-detuned, ^/-polarized beam. 
An alternative scheme would be to use the three beams with parallel polarizations, as 
illustrated in Fig. [T]^. This scheme allows for the use of lower power for the blue-detuned 
beam, by about a factor of three for the parameters that we will consider in section [3} 
but results in larger vector shifts as we will discuss in the next sections. The scheme 
that we propose makes use of four beams with parallel linear input polarizations, as 
shown in Fig. [T]i. The additional blue-detuned beam compensates for the vector shifts 
of its companion blue-detuned beam, as we will show. 

Fig. [2] illustrates a trap generated using the configuration of Fig. and the 
parameters in Ref. p]. About 2000 atoms were trapped in a 1-D lattice with a 50 ms 
lifetime. This advance in the group of A. Rauschenbeutel represents an important 
milestone towards the micro-manipulation of ultra-cold atoms using evanescent field 
traps. 



2.3. HEii mode - electric field polarization 

The fundamental mode HEn is often referred to as "quasi-linear" when excited with a 
linearly polarized input beam. However, for a dielectric waveguide in the strong-guiding 
regime with indices of refraction rii ^ 1.5 inside the waveguide and n2 ~ 1.0 outside, the 
HEii mode actually exhibits a significant ellipticity for a < A/2, leading to vector shifts 
of the Zeeman sublevels. Formally, for a linearly polarized input, the evanescent field 
E = {Ex, Ey, Ez) of the fundamental mode propagating in the fiber can be expressed as 
follows for r > a [Tl lMliS] : 

Ex{r,(j),z,t) =Aiin- — 4a 

2quKi{qua) 

X [(1 - sn)Ko{qnr) cos(<^o) + (1 + Sn)K2{qur) cos(20 - ^oW'-'^'-^'^'K 

Ey{r,(f),z,t) =Aiin- — 4b 

2qiiKi{qua) 

X [(1 - su)Ko{qur) sin(^o) + (1 + Su)K2{qiir) sin(20 - ^o)]e'^'''-^'''\ 
Ez{r, 0, z, t) = ^Au,^i^^iri(gnr) cos(0 - y.o)e^^"*-^"^\ (4c) 

with 
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Figure 1. Trapping schemes, a) Field intensity \E\^ in the plane transverse to 
the fiber for a single, x-polarized beam at A = 937 nm. \E\^ is normalized to the 
intensity just outside the fiber /q = \E(r = — 0)^, with a = 250 nm and 

a+ = a + |r — a| , r — > a. The red dashed arrow indicates the input polarization, b) 
Trapping scheme used in Ref. [B]. Red(blue)-detuned beams are shown by red(blue) 
thick arrows. Input polarizations are shown by the thin arrows. A single, y-polarized 
blue-detuned beam is used, c) Three-beam scheme with parallel x-polarizations. All 
beams have an intensity maximum in the x — y plane along the direction of the input 
polarization, d) A second x-polarized blue-detuned beam is added to compensate for 
the vector shifts, as discussed in the main text. 
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Here, denotes the azimuthal position in the transverse plane (Fig. [T|l), yjg 
indicates the polarization axis for the input polarization relative to the x axis, ni 
and n2 are the indices of refraction inside and outside the waveguide, is the mode 
propagation constant, l/Zin is the characteristic decay length for the guided mode inside 
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Figure 2. A ground state trapping potential J/trap is generated by two orthogonally 
polarized evanescent fields, confining Cs atoms outside a 500 nm diameter optical fiber. 
Input polarizations are denoted by the arrows, a) x-y plane, b) x-z plane, f/trap results 
from two counter-propagating red-detuned beams (1064 nm, 2 x 2.2 mW), and a single 
blue-detuned beam (780 nm, 25 mW), as shown in Fig. ^ [5]. The standing wave 
structure of the attractive red-detuned field and the repulsive force from the blue- 
detuned beam enable 3D confinement of Cs atoms at each minimum of f/trap near the 
dielectric waveguide. J/trap diverges as the surface is approached due to the attractive 
van der Walls force. 



the fiber, l/gn is the characteristic decay length for the guided mode outside the fiber, 
Aiin is the real- valued amplitude for the hnearly polarized input, Ji is the l-th Bessel 
function of the first kind, and Ki is the l-th modified Bessel function of the second kind. 

It is clear from Eq. Q that the electric field intensity is not azimuthally symmetric. 
For a beam polarized along e^,, i.e. (po = 0, the intensity at the fiber's outer surface is 
maximum for = 0, vr and minimum for = ±7r/2. 

Notably, the evanescent modes of the nanofiber have a significant longitudinal 
component along the fiber propagation direction, which is 7r/2 out-of-phase with 
the transverse components {Ex,Ey) (Eq. (4c)). At the outer fiber surface, Ez is 
maximum for = (po,(po + tt (i.e., along the input polarization axis), and vanishes for 
= <y9o±7r/2. For an x-polarized input at 937 nm and a nanofiber of radius a = 250 nm, 
^j^i" (r = a+, = 0) ~ 20%. As a consequence, the polarization of a single propagating 
beam will be elliptical everywhere except for = (/^o ± 7r/2. The ellipticity of the beam 
will be maximum for = y^o, (po + ir as is illustrated in Fig. [3| giving rise to significant 
vector shifts, which we describe in section |3} 

We can re-write Eqs. ^ as follows: 
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Figure 3. Electric field E(a;, of a single propagating beam in the plane y = 0. 

The input beam is x-polarized. The electric field Re[E(a:, j/, z, i)], with E(a;,?;, z,i) 
defined as in Eq. [4] is shown by the blue arrows. The red arrow indicates the beam 
propagation direction. The field is shown for a) = 0, b) ujt = it /2, and c) ut = tt. 



E,(r,(/),z,t) = Ae^('^*-^"^), (6a) 
Eyir, 0, z, t) = ^e^^"'*-^""), (6b) 
E,{r, 0, z, t) = zCe^('^*~'^"^), (6c) 

where A, B, and C are real functions of r and (p. In particular, if one combines a 
forward-propagating beam E*^^™'^) expressed as Eq. ^ with a backward-propagating 
beam of same amplitude and input polarization E'^'^'"'^) = y4e**^'^*+^"^^e2:-|--Be*^'^*+^"^^ey — 
the total field can be expressed as: 

E(tot) _ g(fwd) ^ g(bwd) ^2 [A cos{Puz)e^ + B cos{(3uz)ey + C sin(/3n^)e,] ■ e'^'K (7) 

The resulting electric field E(*°*) = E^^"^*^) +E('^™'^) forms an optical lattice with spatially 
rotating linear polarization as illustrated in Fig. |4j In particular, the polarization state 
of the field rotates between the pure linear x and 2;-polarizations along z at (/) = 0, as 
illustrated in Fig. [5j 

2.4- Cancellation of the vector shifts 

In Ref. |6], the properties shown in Figs. |4] and [s] were used to cancel the vector shifts of 
ground and excited mp states in Cs for the pair of red-detuned trapping beams for the 
configuration in Fig. [ijo. However, vector shifts due to the single blue-detuned beam 
were zero only for = ipo, ipo + n. Although the atoms are trapped at = ipo, ipo + vr, 
inevitable fluctuations of the atom position will lead to non-zero vector light shifts of 
both ground and excited states. 

The scheme in Fig. [T]^ allows for the use of reduced power for the blue-detuned 
beams as compared to Fig. [T)d but with the consequence of large vector shifts from the 
ellipticity of the electric field even for = 0. We will therefore not consider this scheme 
in the next sections. 
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Figure 4. Total electric field E(x,y, for two counter-propagating beams in the 
plane y — 0. The input beams are cc-polarized. The electric field Tie[E{x , y , z , t)] is 
shown by the blue arrows. The red arrows indicate the beams propagation directions. 
The electric field is shown for a) uit = 0, h) ut = tt/A, and c) ujt = tt. As opposed 
to Fig. [3] the polarization of the electric field is linear at any point |r| > a (i.e., the 
polarization vector has no ellipticity and E does not rotate in time at a given position 
r as in [3]). 
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Figure 5. Electric field amplitude after interference E(*°*) = E^^*'^) + E^''^'^) of two 
937 nm beams with 6 ft — 0, Sit t — and r = a+ as in Fig. [TJi (i.e., a:;-polarized inputs 
with (ySo = 0). The fields are normalized to the intensity Iq at r — a^, (j) — 0, z — 0. 
a) axial direction z (at = 0). b), azimuthal direction (at z — 0). In particular, 
j;(tot) Yyg^g a fixed linear polarization at any given point r which rotates as r is varied. 



By contrast, the vector shifts of both the ground and excited states can be canceled 
for both the red and blue-detuned fields by using pairs of counter-propagating beams, 
as shown in Fig. [T]i. In the x — z plane, the vector shift for each pair becomes 

Hi oc (a(i)(u;(^"^)) - a^^\u^^'"'^^))^ with u^^'"'^^ - u^^'"'^\ where u;(f«d,bwd) 

are the 

angular frequencies for the forward and backward propagating beams, and 6 ± Sfb/2 
are their detunings from the atomic transition frequency Ua, with two-photon detuning 
5fb = cj'-^"'^^ — uj^^'^'^K For an atom in the x — z plane, the total electric field is also 
contained in the x — z plane, such that the scalar product ^E'^^^ x E'^+)^ ■ F in Eq. ^ 

is proportional to Fy. 

In the case of the red-detuned lattice, oof^'^'' = ^^el^^ ^^"^ Hi*^^^ = 0, precisely 
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as in Ref. [6]. Adding a blue-detuned lattice with 5fb = would result in two 
superimposed lattices with unmatched spatial periods 27[ //31'f-, 271 //3\]^'^. To avoid this 
effect, the interference between the counter-propagating blue-detuned fields e[,^"^^ and 
-^Wucf '^^^ tie averaged over times short compared to the time scale of the motional and 
internal dynamics of a trapped atom by offsetting the frequencies of the two fields by 
6fb ^ (entrap, ^hfs), whcrc cutrap and (5hfs are the trap angular frequency and the hyperfine 
splitting for the ground state, respectively. This will also suppress spurious two-photon 
processes (e.g., two-photon Stark shift [16]) as well as parametric heating due to intensity 
modulation [17]. 

For w^iue '^"'^'' = + ± Sfb/2), we achieve a vector shift cancellation for the 
blue-detuned field to the first order in 1/6, namely: 

^r'^oc^§ + (9(1/5=^). (8) 
For typical values of 5 = 85 THz and 6fb = 30 GHz, Sfb/S = 3.5 x 10~^. 
2.5. Magic wavelengths for an evanescent field trap 

To make the nanofiber trap state-insensitive, it is necessary to cancel the differential 
scalar shift 5?7scaiar by operating the trap at the magic wavelengths, as proposed in 
Ref. [33], in which only the effects of the scalar and tensor shifts were considered. 
Here we deal with the full complexity of the vector field E(r) and the resulting 
vector light shifts. We numerically determine the red-detuned and blue-detuned 
magic wavelengths, following the procedure described in Refs. [301 EH US, US]. We 
include all the hyperfine levels F and Zeeman sublevels nip of the electronic states 
{6^1/2, ■ ■ ■ , 15^1/2}, {6P1/2, ■ ■ ■ , IIP1/2}, {6P3/2, ■ ■ ■ , IIP3/2}, {61^3/2, ■ ■ ■ , 11^3/2}, and 
{6D5/2, ■ ■ ■ , II-D5/2}. The effect of the tensor shifts on the excited state is manifest in 
the quadratic splitting of the mp' sublevels (Fig. |6|. We find a red-detuned magic- 
wavelength located around 935nm, in accordance with the previously published values 
[311 119] . In the next sections, we will use the value Ared = 937 nm, that cancels 5?7scaiar 
for the 6P3/2 excited state \F' = A^mp' = 0). We choose P' = 4 due to its relevance 
to coherent two-photons processes [151 US [HI- There are several blue-detuned magic 
wavelengths [221112] • For our trap, we use the magic wavelength Abiue at approximately 
687 nm [35] . Since this is the second closest blue-detuned magic wavelength to 852 nm, 
it has the second highest ground-state polarizability and therefore requires the second 
lowest optical intensity to generate the required trapping potential (we do not consider 
the magic wavelength at 792 nm, as it is too close to the 8S'i/2 to 6P3/2 transition at 
794 nm). 

We have neglected higher order processes in our analysis, including two-photon and 
electric quadrupole transitions, near 687 nm [HU] . 
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Figure 6. Magic wavelengths of the Cs D2 line. We display the light shift U\s for a 
linearly polarized beam with constant intensity 2.9 x lO^W/m^ around a), the blue- 
detuned magic wavelength at Abiuc — 687 nm, and b), red-detuned magic wavelength 
at Ared — 937 nm. 



3. Numerical results: trapping potentials 

Using the atomic interaction Hamiltonian in the dipole approximation with the actual 
polarization profile of the evanescent field, we proceed to analyze the adiabatic potentials 
for the nanofiber trap for a Cs atom in its 6S1/2 ground and 6P3/2 excited states. 



3.1. Total potential 

For a specific atomic state of Cs, the total atomic trap potential [/trap consists of the 
total light-shift potential Uis calculated from the full Stark shift Hamiltonian (Eq. ([2])), 
as well as the surface interaction potential of an atom with the dielectric waveguide 
[/surface, namely 



Utra,p Uig -\- [^surface • (9) 

The Casimir-Polder interaction between the atom and dielectric surface has a 
significant effect on the atomic motion at distance scales near 100 nm 1 ^ 1 ^ [HS f [Si] . 
The surface potential of a ground state Cs atom near a planar dielectric surface can be 
reasonably approximated by the van der Waals potential which decays as d~^, where 
d = r — a: 

C3 

[^surface 5 (-^0) 

where we use C^/h = 1.2 kHz /im^ ^T\. Because the retarded Casimir-Polder forces 
decrease faster away from the surface than the van der Waals forces, using f/surface 
overestimates the surface interaction at the trap location d ^ 200 nm. Additionally, 
the curvature of the nanofiber cylindrical geometry reduces the potential strength even 
further [211 ES] . The d~^ scaling of the van der Waals expression for a planar surface 
is therefore an overestimate of the actual surface potential, and is used for simplicity 
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Figure 7. Radial dependence of the trapping potential of the ground and excited 
states for the parameters used in Ref. [6] at z = 0. The polarization configuration 
is the same as Fig. [TJd. The energy sublevels of the ground states F — 3 and F — A 
of 65i/2 are shown as solid green and dashed black curves, and the F' — A sublevels 
of the electronically excited state (6P3/2) are shown as red dashed curves, a) Radial 
potential along (/i = 0. The trap minimum is located about 230 nm from the fiber 
surface. The excited state is un-trapped, and split by the tensor shifts, b) Radial 
potential along (j) = 7r/2. Both ground and excited states are not trapped. The ground 
states exhibit a splitting due to the vector shifts induced by the elliptical polarization 
of the blue-detuned light, and the excited states are shifted by the vector and tensor 
shifts. 



in the calculations presented here, with more complete expressions for Cs presented in 
Ref. [21]. Furthermore, we neglect any dependence on the mp' sublevels of the excited 
state 6P3/2, and simply approximate t/gurfece ^ 2 f/surface [56]. 

We calculate the adiabatic potential of Eq. ^ by diagonalizing the total interaction 
Hamiltonian H = H\s + -/^surface at each point in space, where -f/surface is the scalar 
surface Hamiltonian. At each point r(r, 0, z), we obtain a set of eigenstates and the 
corresponding eigenenergies. These eigenstates are superpositions of the \F,mp) bare 
Zeeman sublevels. Due to the complex polarization of the trapping fields, the energy 
eigenstates are not necessarily eigenstates of any projection of the angular momentum 
operator. 



3.2. Effect of the light shifts in a "non-magic" trap 

First, we consider the trapping parameters for the experiment by Vetsch et al. [0]. 
Despite its impressive experimental success, the ground-state levels exhibit splittings 
that impair the ground state coherence. In the realization of Ref. [6], the two-color 
evanescent trap is constructed using a pair of counter-propagating x-polarized (y^o = 0) 
red-detuned beams E^ed = eJ^^'^^ + eJJIJ'^^ (P^cd = 2 x 2.2 mW) at A^ed = 1064 nm, 
forming an optical lattice, and a single repulsive y-polarized (yjo = ^/2) blue-detuned 
beam Ebiue (Pbiue = 25 mW) at Abiue = 780 nm. The Si02 tapered optical fiber has a 
radius a = 250 nm in the trapping region. 
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Figure 8. Azimuthal dependence of the trapping potential of the ground and excited 
states for the scheme in Fig. and for the parameters used in Fig. [7] r — a = 230 nm, 
z — 0. a) The ground state sphtting is minimum for (/) — and (f> — Everywhere 
else, the polarization of the blue-detuned field induces large vector shifts, b) Expanded 
view of a) near a trap minimum at = tt. The Zeeman-like splitting of the ground 
states is large even for small azimuthal angles. The excited-state level structures are 
greatly altered by the combined vector and tensor shifts. 



Fig. [T] shows the radial trapping potential Utr.a,p{r, (j), z) of the ground states F = 3 
and F = 4 of QS1/2 and excited states F' = 4 of 6P3/2, for 2; = and = (x-axis) (Fig. 
[7^) and for = 7r/2 (?/-axis) (Fig. [TjD). The energy sublevels of the ground states at 
the trap location (0 = 0) are degenerate, as both trapping fields are linearly polarized 
as illustrated in Fig. |4} The excited state energy sublevels are shifted due to the vector 
and tensor shifts. The trap depth for the ground state is t/depth = —0.4 mK, located at 
r — a ^ 230 nm and = 0, whereas the excited states are not trapped at all. 

The azimuthal dependence of the trap potential reveals a significant inhomogeneous 
broadening of the energy sublevels due to the ellipticity of Ebiue for 7^ 0, vr (Fig. [s]). 
To estimate this broadening, we assume that the potential is harmonic around the trap 
minimum. By fitting the ground state F = 3 potential with a harmonic potential 
around = vr, we obtain an azimuthal trapping frequency z^trap^ — 150 kHz. For an 
atom in its azimuthal motional ground state |n)^ = in such a potential, the half- 
width (Tr^ of the corresponding single- atom distribution is given by a^^ = {{^(pY) — 
ftvapCtt, = ^ i^jnut ~ ~ azimuthal half-width of ~ 2°). This leads to 

fast decoherence of the hyperfine and Zeeman levels, even with ground state cooling. 
Specifically, we estimate a spin- wave coherence time = l/6iy^ < 5 /is, derived from 
the Su^ = 200 kHz splitting between the sublevels of the F = 4 atomic ground state 
16 nm away from the trap minimum. This is significantly limited compared to the 
quantum memory performances of atomic ensembles in optical lattices (see, e.g., [57]). 

Finally, we plot in Fig. [9] a cross-section of the axial potentials showing the axial 
confinement of the ground and excited states. 

The excited states are untrapped in all directions except along the fiber axis z for the 
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Figure 9. Axial dependence of the trapping potential for the ground and excited 
states for the scheme in Fig. and for the parameters used in Fig. [7j a) Longitudinal 
potential along (p — 0. The distance from the fiber surface is set to r — a = 230 nm at 
the trap minimum, b) Longitudinal potential along (f) = tt/2. The distance from the 
fiber surface is again set to 230 nm. 



parameters of Ref. [6]. An atom excited to these untrapped potentials will experience 
dipole-force fluctuations, leading to heating [29] and preventing near-resonant driving 
of the optical transition 



3.3. State-insensitive trapping potential 

Now, we analyze our proposed "magic compensation" scheme (as illustrated in Fig. [Ij; 
demonstrating how using magic wavelength beams and compensating the trap ellipticity 
can reduce inhomogeneous broadening of the Zeeman sublevels in a nanofiber trap. For 
this trap, we use a pair of counter-propagating x-polarized {ipo = 0) red-detuned beams 
Ered = E^""^^ + eJ|;J'^^ (Pred = 2 X 0.95 mW) at the "magic" wavelength A^ed = 937 nm, 
forming a 1-D optical lattice. Counter-propagating, x-polarized blue-detuned beams at 
the second "magic" wavelength Abiue = 687 nm are used with a power Pwue = 2 x 16 mW. 
The resulting interference is averaged out by detuning the beams hy Sfb = 30 GHz, as 



explained in section 2.4, leading to a first-order cancellation of vector light shifts as 
expressed by Eq. (|8]). The beam intensities are chosen to generate a trap of similar 
depth as the one demonstrated in Ref. [B]. The resulting adiabatic potential f/trap 
allows for state-insensitive 3D confinement of cold Cs atoms around a Si02 nanofiber of 
radius a = 250 nm. 



In Fig. 10, we show the radial trapping potential f/trap 0, ^) of the ground and 



excited states for z = 0, = (a;-axis) (Fig. 10 1) and for z = 0, = tt/2 (y-axis) 



(Fig. 10 d). Because the trapping fields are now effectively linearly polarized, the ground 
states are degenerate at both = and (p = tt/2. In contrast to a non- magic wavelength 
trap, the excited states are trapped with gradients that closely map that of the ground 
states. The sublevels of QP3/2 are still non-degenerate due to the tensor shifts. For 
-Pred, -Pbiuo Specified above, we find that the trap depth is t/depth = —0.4 mK, located at 



A state-insensitive nanofiber trap 



15 



c) 



0.50 



0.25 



0.00 



-0.25 



-0.50 



0.00 




-0.25 



-0.50 
100 



200 400 600 800 

r - a (nm) 



1 










^=0 



150 



200 250 

r - a (nm) 



300 



1000 



b) 



S 



350 



d) 



0.50 



0.25 



0.00 



-0.25 



-0.50 



0.00 




-0.25 



200 400 600 800 1000 

r - a (nm) 




150 200 250 

r - a (nm) 



350 



Figure 10. Radial dependence of the trapping potentials of the ground and excited 
states using the magic wavelengths and the configuration shown in Fig. [TJi. All beams 
are polarized along ip — (i.e., (po—O). The 937 nm beams each have a power of 0.95 
mW. The 687 nm beams each have a power of 16 mW. a) Radial potentials along <j) — 
(i.e., ifo = 0). The trap minimum for 6S'i/2 is located about 210 nm from the fiber 
surface. Both electronic ground and excited states are trapped, with residual splittings 
of the excited states due to the tensor shifts, b) Radial potential along (j) = 7r/2. c), 
d) Expanded view of a), b) around the trap minimum. 



^trap — a = 210 nm and = 0, vr. 



The azimuthal confinement of tlie atoms is sliown by Fig. 11 In contrast to tlie 
configurations sfiown in Figs. [1]d and [T]^, tfie inliomogeneous Zeeman broadening from 
tlie ellipticity of Ebiue is greatly reduced thanks to the compensation scheme of Fig. [T]i. 
It is non-zero, however, as expressed by Eq. ([s]). The remaining sphtting of the F = A 
ground state is 5^ ~ 700 Hz, limiting the coherence time to r < l/(5jy = 1.4 ms. 

In the case of perfect cancellation of the vector shift with = 0, a residual 
non-zero ground state splitting dv^f, would still arise from the different scalar dynamic 
polarizabilities of the 65*1/2 F = ?> and F = 4 ground states [21]. For atoms in their 
azimuthal motional ground state \n)^ = jO),/,, the single-atom distribution half-width is 
— 30 nm (or ~ 4°) with azimuthal trap frequency t'trap^ — 44 kHz obtained from 
a harmonic fit of the potential around (p = it. We estimate a spin- wave coherence 
time Tm = 1/A((5z/^) < 30 ms, coming from the spread A (^z/^) = 5//^ (0 = vr) — 
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Figure 11. a) Azimuthal dependence of the trapping potential of the ground and 
excited state for the "magic compensated" trap (Fig. [T]i) with the parameters used 
in Fig. 10 for r — a = 210 nm. b) Expanded view of a) around a trap minimum. 



The effect of the compensation beam in Fig. [Tp is to suppress the vector shifts and to 
reduce the ground-state splittings Siy^ in the transverse plane for </> 0, tt. 



6^^ (0 = TT + (T(^) ~ 30 Hz of the atomic ground states for the F = 3 — )■ F = 4 transition 
frequency. 

We note that the longest achievable coherence time in the "magic compensated" 
adiabatic potential in the absence of ground-state splitting 5//,^ would be limited by 
spontaneous Raman scattering driven by the trapping beams 
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Figure 12. Axial dependence of the trapping potentials of the ground and excited 
states for the "magic compensated" trap (Fig. [TJl) with the parameters used in Fig. 10 
a) Longitudinal potential along (/) — for r — a = 210 nm. b) Longitudinal potential 
along (f) = 1^/2 for r — a — 210 nm. 



Finally, we also plot the axial potentials in Fig. [T2| showing the confinement for 
both the ground and excited states. 

Thanks to the use of the magic wavelengths, the excited states are trapped in 
all directions. This results in greatly suppressed dipole-force fluctuations, allowing for 
on-resonance driving of the optical transition. 
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4. Conclusion 

We have proposed and analyzed in detail a scheme for a state-insensitive optical 
nanofiber trap which utilizes realistic experimental parameters. The "magic" trapping 
wavelengths of 937 nm and 687 nm for Cs atoms are readily available using 
semiconductor laser sources, and require a reasonable power for trapping. Tapered 
optical fibers can be made with sub- wavelength diameters and high transmission, as 
has been demonstrated experimentally [601 [61]. In R-sf. [62] we describe an experiment 
for trapping Cs using Fig. [T]i configuration and explicitly demonstrate features of the 
"compensated magic" trap for a nanofiber. 

Furthermore, extension to other nanoscopic dielectric waveguides [63] would make 
evanescent optical trapping possible in integrated hybrid quantum devices [M]. It is 
worthwhile to note that the compensation scheme of the vector shift would work at any 
wavelength, increasing the ground state coherence time in a straightforward manner. 

The proposed "magic compensated" scheme allows for in-trap resonant processes, 
and leads to increased ground state coherence time relative to uncompensated schemes. 
These properties will make quantum-state engineering more feasible in such a trap, 
allowing for a wide range of experiments including creating quantum memories, coupling 
of single atoms and ensembles to optical or mechanical resonators, and studying 1-D 
spin chains. 
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